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Abstract 

The hierarchy of integrable equations are considered. The dynamical 
approach to theory of nonlinear waves is proposed. The special solutions 
(nonlinear waves) of considered equations are derived. We use powerful 
methods of computer algebra such as differential resultant and others. 

1 Introduction 

The paper is organized as follows: section 1 describes the basic facts of hierarchy 
of nonlinear equations, Baker-Akhiezer (BA) function theory for finite-gap in- 
tegrationn method. Particular solutions in terms of generalized Novikov, Lame, 
Hermite type polynomials are presented. 



2 Preliminaries 

We review in this section some basic facts about Baker- Akhieser function which 
will be used in the sequel. For systematic treatments and proofs, we refer the 
reader to §, ||. 



[9,-C/,9t-yW] =0, (2.1) 
U = U{x,t,X), y(^) =y(^)(a;,i,A), (2.2) 

where V'^^^ — VqX^ + ViA*^^^^-* + . . . Vn- The nonlinear equations in the form 
(2.2) are obtain as compatibility condition of the following linear differential 
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equations 



5,* = C/(a;,i,A)*, 



(2.3) 



where 4' is a vector-collumn (^i . . . 5'„)"^. The system (2^) play an impotran 
role in the theory of finite-gap integration method. The function which is 
general solution of the system (2.3), is meromorphic function on the Riemann 
surface and have esential singularities of prescribed form near fixed points on 
Riemann surface. This surface is defined by 

det(i(x,i,A) -/il) =0. (2.4) 

We call finite-gap solutions of ( |2.2| ) these for which there exists the meromorphic 
in A matrix function L(x, t, A) such that 



= [C/,L], L 



Lt^[V\L], 
= ioA^+LiA^-i 



(2.5) 
(2.6) 



3 Resultant of two differential operators 

Let us consider the following two differential operators 

n m 

A = Y,akD\ a„7^0, B = Y,biD\ b„, ^ 0, ak^keCf^ (3.7) 

fc=0 (=0 

where D = d/dx, Cj*-'-* is the set of fc(/)-th differentiable functions on the 
interval / of real variable x. Right resultant of operators A and B RRes {A, B) 
is called the determinant of resultant matrix R of degree (m -I- n) . Analogically 
left resultant of operators A and B LRes {A, B) is called the resultant matrix 
R* of degree {m + n), i.e. RRes (A, B) = det(i?), LRes (A, S) = det(i?*). By 
definition LRes (A, S) = RRes (A*, B*), where A* and B* are the conjugated 
operators. The differential resultant answers to the question when the operators 
A and B have right (left) divisor i.e. when the overdetermined system Ay — 0, 
By — (or A*y = 0, B*y — 0) has soluion. By construction to define right 
resultant RRes {A, B) we act with operators Z)™~^, . . . ,D, — 1, by left to A 
and with operators D"~^, _D" ~ 1. to B. In this case the system Ay = 0, 

By = has the following form 

n+s m+p 

$]aM2/^'^=0, s = 0-(m-l), ^ 6,,^/) = 0, p = - (n - 1),(3.8) 

fe=0 1=0 

where the coefficients au.s and hi^p are computed by: 

«M-E(D«-'' ^'^^ = E(-)fc^'- (3-9) 

i=0 3=0 ^ ^ 
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The system (3.8) is constructed by homogeneous hnear algebraic equations in 
terms of the foUowing variables ■ ■ ■ ,y' , y- The matrix of coefficients 

of this system is the right resultant of the matrix R, det(i?) = RRcs {A, B). 



RRes (A, B) 



'^n+^Ti-~l,m — 1 ^n+m— 2,m— 1 
^n+m— l.m— 2 










^n+m — 2,71—1 
bn+m-2,n-2 



&0,Ti-l 

6o,o 



.(3.10) 



If the system Ay = 0, By = is consistent then dct(i?) = RRes (A, B) = 0, in 
the opposite case the system is not consistent. 

For applications of differential resultant see for example [?] . The most im- 
portant of them are the following: the cirterion of existing the greater right 
and left divisor of the system of operators; the criteion of concictency of the 
system of linear ordinary differential equations with one unknown function; the 
criterion of commutation of two linear differential operators; he criterion of ex- 
isting polynomial and exponential solutions of linear differential equtions with 
variable coefficients; the criterion of factrorization of operators of degree n into 
products of operators of degrees n — 1 and 1. 



Example. Let us introduce the following equations, 

Ll = Ay^ x'^y" + xy' - {x'^ + l/4)y = 0, 
L2 = By^ 2xy" + (3 - 4x)y' + {2x - 3)y = 0. 



(3.11) 
(3.12) 



Using ( p.lO| ) we obtain 



x^ 


3x 


-a;2 +3/4 


-2a; 







X 


-a;2 - 1/4 


2x 


5a; - 4 


2a; - 7 
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2a; 


3 -4a; 


2x - 3 



RRes (A, B) 



Problem 1. Let us define two LODE: 

n m 

Ay = ^ flfei^^j/ = 0, a„^0, By^J2biD'y^0, ^ Ca^, 6; G ^^.13) 

k=0 1=0 

Find differential resultant of operators A and B. 

Algorithm 1. We solve problem 1 by the procedure DIFRESULT(a, n, b, m, x); 
Input 
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a is the array of coefficients of differential operator A] 

n is the degree of differential operator A; 

b is the array of coefficients of differential operator B; 

m is the degree of differential operator B; 

X is independent variable. 

Output: 



Using (3.£) with given coefficients of differential operator A and B we com- 
pute the ellements of the resultant matrix , after that we we find the determi- 
nant. The output is RRes [A, B) = det(i?), where R is the resultant matrix. 

4 The generalized Riccati equation 

4.2.1. Prime right divisor of operator L. 

Let us define the differential operator L of degree n: 

n 
k=0 

A necessary and sufficient condition for the following factorization: 

L = L2{D — a), ord_L2 = n — 1 
is given by RRes (L, D — a) = 0. 



RRes [L, D-a) 



1 ~a 






02 

71 — 3 

1 





ai 



2) 



-a 
1 



oo 



-a 
-a 



Then we find the following generalized Riccati equation of first kind N{a) of 
degree (n — 1): 

N{a) = L{a) + M(a), 

where 

n 

.k 



L{a) = ^Ofca* 

fc=0 

and L{a) is called he pseudo-characteristic Riccati equation of first kind 

n 

M(a) = ^flfeMfe^i, Mk-i = {D + a) 



Q(.4.14) 



fe=0 
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and M{a) is called the reduced Riccati equation of first kind. When n = 2 the 
generalized Riccati equation of first kind have the form 

A^(q;) = a2a' + (1,20^ + aia + ao 

4.2.2. Prime left divisor of operator L. Let us define the following 
differential operator L of degree n: 

11 

L = Y^ akD'', a„ 7^ 0, afe G 

k=0 

To find factorization of following type, 

L = L2{D — a), ordI/2 = n—1 

the necessary and sufficient condition is RRes {L, £) — a) = 0, or to find factor- 
ization of the following type, 

L = [D — a)Li, ord Li = n — 1 

the necessary and sufficient condition is LRes (L, D — a) = RRes {L*, D + a), 
i.e. 



a 



Jn-2) (r, 



a 
1 



a 
a 



Thus we obtain the generalzed Ricati equation N*{a) of second kind and degree 
(n-1): 

N*{a) = L*{a) + M*{a), 

where 

L*ia) ^rL = Y: E(-l)' , = ^^(^ + (4-15) 

r=0 k=r ^ ^ 

Here r is the conjugation operator, L*{a) is pseudo-characteristic Riccati equa- 
tion of second kind and 

M*{a) = E ^(-1)'= 4'=-'-^M,_i, Mr-, = {D + a)('-i) - a^'-i^^lS) 

r=lfe=r ^''^ 

M* (a) is the reduced Riccati equation of second kind. 

When n = 2 the generalized Riccati equation of second kind have the form, 
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N{a) = a20e' + a^o? + (— ai + 2(i.^q. + ao — + = 
As a consequence of problem 1 we have: 

Problem 2. Let us define the following equations, 



n 



Lyy = ^ akD'^y = 0, o„ ^ 0, ak € 



(4.17) 



fc=0 



r=0 k=r ^ ^ 



t'\' = 0. 



(4.18) 



Find the generalized equation of first kind, 



Ni{a) =RRes(Li,D~a) 
and the generalized Riccati equation of second kind, 

AT* (a) = RRes(L^,L» + a)). 



Algorithm 2. We solve the problem 2 by the procedure 
RICCATI{a,n,pp,x). 

Input 

a is the array of coefficients of differential operator L\; 

n is the degree of differential operator Li ; 

pp is the kind of seeking Riccati equation [pp = 1 and pp = 2); 

X is the independent variable. 

Output: The generalized Riccati equation of first kind {pp = 1) and of second 
kind [pp = 2). 

Rl: If pp = 1 then go to R2 else go to R3; 

R2: <C i2 = (£) - a)y; 

b is the array of coefficients of operator L2; 

m := 1 ; Ni{a) := DIFRESULT(a, n, b, m, x); return Ni{a) >; 
R3: <C = tLi, t is operator of conjugation. 



a* is the array of coefficients of differential operator if, L2 = {D + a)y; 

b* is the array of the coefficients of differential operator Lji 

m := 1 ; N2{a) := DIFRESULT(a*, n, b*, m, x); return N2{a) »; 
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5 ODE resolved by algebraic means 



4.3.1. Linear differential equations with exponential solutions of type 

y = exp(— ax), a = const. 

The necessary and sufficient condition for existence of sucli solutions is the 
pseudo-characteristic Riccati equation of first kind L{a) and of second kind 
L*{a) to have solution a = const, i.e. RRes [L.D + a) =0 or RRes (L*, D — 
a) = 0. Then the linear differential equation Ly = X]fc=o o-kD'^y have the 
form Ly = L2{D + a)y, where L2 is differential operator, ord(Z/2) = n — 1, 
and the differential equation have the exponential solution of following type 
y = exp(— aa;). In the case of linear equation of second degree 

Ly = a2{x)y" + ai{x)y' + ao{x)y = 0, Ofe ^ 0, UkeCj, k = 0^2, 

the problem is to find the coefficients of factrization ai, in Ly = {D — ai){D — 
Oi2)y = 0, where ai = const and a2 = const, i.e. we consider the classes of 
equations for which the associated Riccati equation is of first type or of second 
type in terms of a have constant solutions. Problem 3. Given LEDE of degree 

n, 

n 

Ly = J2akD''y = 0, a„^0,OfeeCf 
fe=o 

Find exponential solution y = exp^ —ax), where a = const.. 

Algorithm 3. In our program the problem 3 is solved by the procedure 
LIVDIF{a,d, n, x). 
Input 

a is the array of cocfhccnts of operator L; 
d is nonhomogeneous part of equation Ly; 
n is the degree of differential operator L; 
x is the independent variable. 

Output If wc find the particular solution of the following form y = cxp(— ax), 
a = const and differential equation is of degree two the coefficients of 
factorization ai, a2 and the fundamental system of solutions yi{x), y2{x) 
arc obtained. If the degree of LODE is n > 2, then the degree of LODE 
is reduced by 1. If the exponential solution does not exists the message is 
received. 

Dl: Algorithm 2 for finding the generalized Riccati equation of first type Ni (a), 
iVi(a) := RICCATI{a, n, 1, x); 
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D2: Li{a) — X]fe=o '^kd^] Li{a) - pseudo-characteristic equation of first kind; 

D3: Ml (a) = Nx{a) — Li(a); Mi{a) - reduced equation of first kind; 

D4: Subalgorithm for finding constant solution of algebraic 
equation Li{a) in terms of a. 
Output: constant solution ao, or 

obtain message that the solution of this type does not exists, 

D5: If ao exists and is solution of Mi{a) — then go to Dll else go to D6; 

D6: Algorithm 2 for finding the generalized Riccati equation of second type 
N2{a). 

N2{a) := RICCATI(a, n, 2, x); 

D7: 

r=0 fe=r ^ ^ 

(L2(Q!)-pseudo-characteristic equation of second type) 

D8: 

M2{a) = N2{a) - L2{a) (5.20) 
(M2(a)-reduced Riccati equation of second type) 

D9: Subalgorithm of finding constant solution of algebraic equation L2{a) in 
terms of a (Output- const, the solution ao, or the message that such a 
solution does not exists) 

DIO: If ao exists and is solution of M2(a) = then go to Dll else go to D13; 

Dll: If n = 2 then go to D12 else go to D13; 

D12: ii a2 = ao; al = -(ai -|- 02); 

yi = exp{j aidx), y2 = yi j exp{j {a2 - ai)dx)dx; (5.21) 

Message "The coefficients of factorization ai and a2 and fundamental 
solutions yi and 2/2 are obtained " H; 

D13: ii yi = exp(/ a^dx), Message "The degree of differential equation is re- 
duced by 1"; New Lode := Sub(?/ = yi J zdx,Ly); In equation 

n 

Ly = Y,akD''y = Q, a„ 0, e (5.22) 

fe=0 

we make the change of variables y = yi J zdx where z is new variable) 
Return New Lode a; 
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D14: Message "The solution of exponential type does not exists" ; Subalgorithm 
for finding of constant solution of algebraic equation L{a) in terms of a. In 
our program this subalgorithm is realized by procedure HAREQ(L(q!), a,x); 
Output: L{a) is algebraic equation in term of a. x is independen variable. 
Output: We have const, is the solution ao of equation L{a), or message, 
such a solution does not exists Example: 

Ly = y" + (8 + sm{xf)y' + 8 sin(a;)^y = 0; (') = d/dx (5.23) 

In the array a(n), n = 2, we write the coefHcients of the equation 

a(2) := 1; a{l) := S + iim{xf ; a(0) := 8 sin(a;)^; n := 2; d := 0; 

(5.24) 

Using the procedure LIVDIF(a, d, n, x); we seek particular solution of the 
following type y = exp^ —ax), where a = const., and if the differential 
equation is of degree n = 2 , then we find the coefBcents of factorization 
ai, a2 and the fundamental system of solutions yi{x), y2{x), also the 
general solution. Then we have the solution of differential equation of 
type y = exp^ —8a;), and the coefRcents of factorization are 

ai =8; a2 = sin(a;)^, (5.25) 

and due to the fact that the equation is of second degree we may found 
and second solution 

2/2 = 2/1 j exp( j{a2 - ai)dx)dx, (5.26) 
and general solution 

y = Clyi+C2y2. (5.27) 

4.3.2. LODE with polynomial solutions To exists in the set of solu- 
tions of the equation 

n 

Ly = J2akD''y = 0, a„ 0, ak G (5.28) 

fc=0 

polynomial of degree p ^ m the necessary and sufficent condition is 
RRes(L, Li™+i) = 0, where 

n 

L = Y^ akD^, D = d/dx, = dT+^/dx, (5.29) 

fc=0 

i.e. using the differential resultant of operators L and we find the 

degree of existing polynomial solution of LODE and after that by metod 
of indefinite coefficents we find the coefficents of this polynomial and if the 
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degree of LODE is n = 2 we may obtain the general solution and when 
n > 2 we reduce by 1 the degree of LODE. 

Problem 4. Given LODE of degree n: 

n 

Ly = J2akD''y = 0, a„ ^ 0, ak € (5.30) 

k=0 

find the polynomial solution. 

Algorithm 4. To solve the problem 4 we apply the procedure 
POLDIF(a,n,x,hs); 

Input 

a is the array of coefficents of LODE Ly = 0, 
n is the degree of differential operator L, 
X is independent variable, 

hs is the maximal degree of the polynomial solution, 
Output: 

If there exists the polynomial solution of degree p < hs, then in the case of 
n = 2 we find the general solution, and when n > 2 the degree is reduced by 
1. If there is no solution we have the message that polynomial soluion does not 
exists. 

PI: i:=l; 

P2: Algorithm 1 for finding the differential resultant RRes (L, -D'); 
P3: If RRes (L, D') = 0, then go to P6 else go to P4; 
P4: If i > hs then go to P12 else go to P5; 
P5: < i := i + 1; go to P2 >; 
P6: P:= i - 1; 

P7: Message "P is the degree of polynomial solution"; 

P8: Subalgorithm for finding the coefficents of the polynomial solution ( Output- 
polynomial soluition) 

P9: If n = 2 then go to PIO else go to Pll; 
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PIO: Find the fundamental system of LODE. 
<C yi = plynomial solution; 



2/2=2/1 J exp ^- J aidx^ ^dx 



Return ?/2,2/i »; 

Pll: The degree of LODE is reduced by 1. 

\l Ui = plynomial solution; 

New Lode := Sub(y = yi J zdx, Ly), where in equation 

n 

Ly = Y,akD^y = ^, «n 7^ 0, ak^C^i (5.31) 

fe=0 

we make the change of variables y = yi J zdx where z is the new variable 

Return New Lode H; 

P12: Message "The finded polynomial solution of the LODE "; 

Subalgorithm for finding the coefficents of the polynomial solution of the 
LODE. 

In our program this subalgorithm is realized by the procedure DEPOLDIF{a, n,x,p); 
Input: 

a is the array of coefBcents of of the differential equation y = 
n is the degree of differential operator L, 
X is independent variable, 

p is the degree of polynomila solution oi Ly = 0. 

Output: The coefficents of the polynomial solution of by the method of indef- 
inite coefficents. 

Example. 

Ly = x'^{\n{x) - l)y" -xy' + y^0, (') d/dx 
The array a(n), n = -j- 2, contain the coefficents of the given LODE. 

a(2) := x^(ln(x) - 1); a{l) := -x; a(0) := 1; n := 2; hs := 3; 

By using the procedure POLDIF{a, n, x, hs) we find the fundamental system 
of solution yl := x;y2 := — ln(x). Then the general solution is given by: 
y :=Cla; + C21n(a;); 
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4.3.3. The equations in terms of exact differentials. Let us define 

n 

Ly = Y^ a^D^y = 0, a„ 7^ 0, afe e C^, D = d/dx (5.32) 



fe=0 

This equation is expressed in terms of exact differentials of there exists the 
following form 

n-l 

Ly = DY^ hkD^y = 0, D = d/dx, (5.33) 

fe=0 

where 

an = bn-i{x), ao = 6o(a;), au+i = {bk+h'k+i{x)) , {') = d/dx, k = 0^{n- 

(5.34) 

The necessary and sufficient condition the equation = to be in exact 
differential is, 

n 

^(-l)'=4(x) = 0, (5.35) 

fe=0 

If the equation 



Ly = J2akD''y = 0, (5.36) 



k=0 



is expressed in terms of exact differentials in the case n = 2 we may find the 
general solution of this equation and when n > 2 we may find fist integral. 
Problem 5. Given LODE of degree n: 



Ly = Y,(^kD''y = fix), a„ 9^ 0, ak G (5.37) 

fe=0 

Check is this solution of LODE in exact differentials and if the equation is of 
this type find the general solution (when n = 2) and first integral (when n > 2). 

Algorithm 5. To solve the problem 5 we apply the procedure 

EQDIF(a,n,dl,x); 

Input: 

a is the array of coefficents of LODE 

n is the degree of differential operator L, 

dl is the nonhomogeneous part of given LODE 

x is independent variable, 

Output: 
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Check is given equation in exact differentials and if it is true find the general 
solution in the case n = 1 and first integral if n > 2, and if it not the case, 
message that the equation is not in exact differentials 

El: 

n 

D = 0, (5.38) 

fe=0 

E2: If D = then go to E4 else go to E3; 

E3: Message: "This equation is not in exact differentials". 

E4: ii h{n - 1) a(n); For k -.^ {n - 2) Step (-1) Until Do 

6(fc) :=a(fc + l)-6(/c + l)', {l)=dldx (5.39) 

E5: If n = 2 Then Go To E6 Else Go To E7; 
E6: 



y^cxp{-J a{l)dx){c2 + J {J dldx + Ci)exp{J a{l)dx)dx); (5.40) 
( Ci , C2 - are constant of integration, y - general solution of LODE) 

E7: 

n-l 

L^y ^ DY.bkD'^y, D = d/dx, (5.41) 

Appendix 3: Differential resultant and algebra of commuting diffren- 
tial operators. 

Example 1 



H^l{pl+pl) + lqiql+ql (5.42) 



We associate to this Hamilonian ( 5.42| ) the following linear differential operators 



L = D^+qi, (5.43) 

M = 16D^ + AOqiD^ + eOpiD^ - {UOql + 2bql)D 

-60(?ipi - 15g2P2, (5.44) 

This example show the abilities of the differntial resultant to find the first in- 
etgrals and algbraic curve associated to given linear diferential operators. Using 
REDUCE 3.4 sintactis we have 
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depend qql , x ; depend qq2 , x ; depend pp 1 , x ; depend pp2 , x ; 

a(0) : =-60*qql*ppl-15*qq2*pp2-zz ; 

a(l) :=-120*qql**2-25*qq2**2; a(2) :=60*ppl; 

a(3) :=40*qql; a(4):=0; a(5):=16; 

b(0) :=qql-ll; b(l):=0; b(2):=l; 

b(3) :=0; 

out nik52kdv; 

res52 : =dif result (a , 5 , b , 2 , x) $ 

off mcd; res52:=num(res52)/den(res52) ; 

write "denumerator:=" ,den(res52) ; 

pl:=deg(res52,zz) ; write "pl:=",pl; 

p2:=deg(res52,ll) ; write "p2:=",p2; 

for i:=0:pl do 

« 

dl(i) :=coeffn(res52,zz,i) ; 

LET df (qql , x , 2) =-3*qql**2-qq2**2/2 ; 

LET df (qq2,x,2)=-qql*qq2; 

LET df (ppl , x) =-3*qql**2-qq2**2/2 ; 

LET df (pp2,x)=-qql*qq2; 

LET df (qq2,x,3)=df (-qql*qq2,x) ; 

LET df (qq2,x)=pp2,df (qql,x)=ppl; 

LET df (qq2,x,5)=df (-qql*qq2,x,3) ; 

LET df (ppl,x,2)=df (-3*qql**2-qq2**2/2,x) ; 

LET df (pp2,x,2)=df (-qql*qq2,x) ; 

LET df (pp2,x,4)=df (-qql*qq2,x,3) ; 

Determineint of resulteint matrix 

i-degree of zz:=0 

j - degree of 11:=0 

d2(0,0) :=0 

i-degree of zz:=0 

j - degree of 11:=1 

2 2 2 4 

d2(0,l) :=4*qql *qq2 - 8*qql*pp2 + qq2 + 8*qq2*ppl*pp2 
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i-degree of zz:=0 

j - degree of 11: =2 

3 2 2 2 

d2(0,2):= - 32*qql - 16*qql*qq2 - 16*ppl - 16*pp2 

i-degree of zz:=0 
j - degree of 11: =3 
d2(0,3) :=0 

i-degree of zz:=0 
j - degree of 11: =4 
d2(0,4) :=0 

i-degree of zz:=0 

j - degree of 11: =5 
d2(0,5) :=-256 

p3:=0 

i-degree of zz:=l 
j - degree of 11:=0 
d2(l,0) :=0 

p3:=0 

15 



i-degree of ■z.-z.:=1 



j - degree of 11:=0 
d2(2,0) :=1 

The associated alegebraic curve has the form (the results are obtaine dby 
differential resultant 

= -256A^ - ^2E)? + %K\ (5.45) 

where 

K = q2PiP2 - qipl + qt + l/^qjqi 

Example 2 (Fordy example [?] 

a(0):=-zz; a(l) : =-5*qq2**2 ; a(2) : =15*ppl ; 
a(3) :=15*qql; a(4):=0; a(5):=9; 
b(0):=-ll; b(l):=qql; b(2):=0; 
b(3):=l; 
out nik53sk; 

res53 : =dif result (a , 5 , b , 3 , x) $ 
off mcd; res53:=num(res53)/den(res53) ; write "denumerator :=" ,den(res53) ; 
pi :=deg(res53,zz) ; write "pl:=",pl; 
p2:=deg(res53,ll) ; write "p2:=",p2; 
for i:=0:pl do 
« 

dl(i) :=coeffn(res53,zz,i) ; 

LET df (qql , x , 2) =-qql**2/2-qq2**2/2 ; 

LET df (qq2,x,2)=-qql*qq2; 

LET df (ppl , x) =-qql**2/2-qq2**2/2 ; 

LET df (pp2,x)=-qql*qq2; 

LET df (qq2,x,3)=df (-qql*qq2,x) ; 

LET df (qq2,x)=pp2,df (qql,x)=ppl; 

LET df (qq2,x,5)=df (-qql*qq2,x,3) ; 

LET df (ppl,x,2)=df (-qql**2/2-qq2**2/2,x) ; 

LET df (pp2,x,2)=df (-qql*qq2,x) ; 

LET df (pp2,x,4)=df (-qql*qq2,x,3) ; 

dl(i):=dl(i); 



(5.46) 
(5.47) 
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•/.write "dl(",i,"): = ",dl(i); 

p3:=deg(dl(i) ,11) ; write "p3:=",p3; 
for j :=0:p3 do 

« d2(i, j) :=coeffn(dl(i) ,11, j) ; 

write "i-degree of zz:=",i; 

write "j - degree of ll:=",j; 

write "d2(",i,",",j,"):=",d2(i,j); 

write >> >> ; 



Determincint of resulteint matrix 
denumerator : =1 



i-degree of zz:=0 
j - degree of 11 : =0 
d2(0,0) :=0 

i-degree of zz:=0 
j - degree of 11:=1 

4 2 2 4 2 

d2(0,l):= - 9/4*qql *qq2 - 3/2*qql *qq2 - 9*qql *qq2*ppl*pp2 

6 3 2 2 

- l/4*qq2 - 3*qq2 *ppl*pp2 - 9*ppl *pp2 

i-degree of zz:=0 
j - degree of 11: =2 
d2(0,2) :=0 
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i-degree of zz:=0 

j - degree of 11: =3 

3 2 2 2 

d2(0,3) :=27*qql + 81*qql*qq2 + 81*ppl + 81*pp2 

i-degree of zz:=0 
j - degree of 11: =4 
d2(0,4) :=0 

i-degree of zz:=0 
j - degree of 11: =5 
d2(0,5) :=-729 

p3:=0 

i-degree of zz:=l 
j - degree of 11:=0 
d2(l,0) :=0 

p3:=0 

i-degree of zz:=2 
j - degree of 11:=0 
d2(2,0) :=0 
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p3:=0 



i-degree of zz:=3 
j - degree of 11:=0 
d2(3,0) :=1 

The associated algebraic curve has the form (the results are obtained by differ- 
ntial resultant) 

= 729A^ - IQ2E\^ + K^X, (5.48) 

where 

E=\{p\+pl) + \ql^ + \qiql (5.49) 
3 1 

K = 3pip2 + ■^q2ql + 2«29i, (5.50) 

printC'A program for calculation differential resultants'): 
print ('written by N.A. Kostov, Z.T. Kostova 8 August 1994'): 

#find the elements of series 
procof f :=proc(s ,k,m,b,x) 

local i,sl,s2,cr; 
si :=max(0,k-m) ; s2 : =min(s ,k) ; 
if k < or k > (m+s) then cr:=0 else 

cr:=0; for i from si to s2 do 
if (s-i)=0 then cr :=cr+binomial(s,i)*b(k-i) else 
cr : =cr+binomial (s , i) *dif f (b (k-i) , x$s-i) ; 

fi; 

od; 

fi; 

RETURN (eval(cr)) 
end: 

dif result : =proc(a,n,b,m,x) 

local il,i2,i2, j2,i3, j3,k,s,l,p; 

#** Prograimn for finding ** 
#** differential resulteint ** 
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r : =array ( 1 . . n+m , 1 . . n+m) : 
for s from to (m-1) do 

for k from to (n+s) do 
ak(k,s) :=procoff (s,k,n,a,x) ; 
od; 
od; 

for p from to (n-1) do 
for 1 from to (m+p) do 

bk(l,p) :=procoff (p,l,m,b,x) ; 
od; 

od; 

for il from 1 to m do 

for jl from 1 to (n+m) do 

if il > 1 and jl < il then r[il,jl]:=0 else 
r[il, jl] : =ak(n+m-j 1 ,m-il) ; 

fi; 

od; 

od; 

for ±2 from 1 to n do 

for j2 from 1 to (n+m) do 

if i2>l and j2<i2 then r [m+i2, j2] :=0 else 
r[m+i2,j2] :=bk(n+m-j2,n-i2) ; 

fi; 

od; 

od; 

print ('"Element of resultant matrix'"); 

# for 13 from 1 to (n+m) do 

# for j3 from 1 to (n+m) do 

# print('r:=' ,r[i3,j3] ,13, jS) ; 

# od; 

# od; 

with(linalg,det) ; 
rresl : =det (r) : 
RETURN(rresl) 
end: 

#alias ( dql=dql (x) , ql=ql (x) , q2=q2 (x) , dq2=dq2 (x) ) : 
#b(0) :=dql-z; 

#b(l):=2*ql; b(2) :=0; b(3) :=1 ; 

#a(0) :=-5*q2*dq2-60*ql*dql-ll; 

#a(l) :=-120*ql**2-(35/4)*q2**2; a(2) :=45*dql; 

#a(3) :=30*ql; a(4):=0; a(5):=9; 

#res3:=difresult(a,5,b,3,x) ; 
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#alias(UU=uu(x) ,WW=ww(x)) : 

#alias (UUx=dif f (uu(x) ,x) ,UUxx=dif f (uu(x) ,x,x) ,UUxxx=dif f (uu(x) ,x,x,x) ) : 

#alias(WWx=dif f (ww(x) ,x) ,WWxx=dif f (ww(x) ,x,x) ,WWxxx=diff (ww(x) ,x,x,x)) : 

#alias(WWxxxx=diff (ww(x) ,x,x,x,x) ,WWxxxxx=dif f (ww(x) ,x,x,x,x,x)) : 

#b(0) :=-Z+diff (ww(x) ,x) ; b(l):=2*WW; b(2):=0; b(3):=l; 

#a(0) :=20*WW*diff (ww(x) ,x)+10*dif f (ww(x) ,x$3)-L; 

#a(l) :=20*WW**2+35*diff (ww(x) ,x$2) ; a(2) :=45*dif f (ww(x) ,x) ; 

#a(3) :=30*WW; a(4):=0; a(5):=9; 

#difresult(a,5,b,3,x) ; 

aliasC DC!l=dql(x) ,Ql=ql(x) ,Q2=q2(x) ,DQ2=dq2(x) ) : 

alias (DQlx=diff(dql(x) ,x) ,DQlxx=dif f (dql(x) ,x,x) ,DQlxxx=dif f (dql(x) ,x,x,x)) : 
alias(DQ2x=dif f (dq2(x) ,x) ,DQ2xx=dif f (dq2(x) ,x,x) ,DQ2xxx=dif f (dq2 (x) ,x,x,x)) : 
alias (DQlxxxx=dif f (dql (x) ,x,x,x,x) ,DQlxxxxx=dif f (dql (x) ,x,x,x,x,x)) : 
alias (DQ2xxxx=dif f (dq2(x) ,x,x,x,x) ,DQ2xxxxx=diff (dq2(x) ,x,x,x,x,x)): 

alias(Qlx=dif f (ql(x) ,x) ,Qlxx=diff (ql(x) ,x,x) ,Qlxxx=diff (ql(x) ,x,x,x)) : 
alias (Q2x=diff(q2(x) ,x) ,Q2xx=diff (q2(x) ,x,x) ,Q2xxx=dif f (dq2(x) ,x,x,x)) : 
alias (Qlxxxx=diff(ql(x) ,x,x,x,x) ,Qlxxxxx=diff (ql(x) ,x,x,x,x,x)) : 
alias (Q2xxxx=diff (q2(x) ,x,x,x,x) ,q2xxxxx=dif f (dq2(x) ,x,x,x,x,x)) : 
#alias (Q2Nxxx=dif f (Q2xx (x) , x) ) ; 
#alias(C12xx=C12xx(x)) ; 
DQl:=Qlx: 
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DQ2:=Q2x: 



b(0) :=DQ1-Z; 

b(l):=2*Ql; b(2) :=0; b(3) :=1; 

a(0) :=-5*Q2*DQ2-60*Ql*DQl-L; 

a(l) :=-120*Ql**2-(35/4)*Q2**2; a(2) :=45*DQ1; 

a(3):=30*Ql; a(4):=0; a(5):=9; 

ppp:=difresult(a,5,b,3,x) : 

rrr :=expand(subs(DQl=Qlx,Qlxx=-4*Ql~2-Q2'2/4,Q2Nxxx=Q2xxx, 
Q2xx=-Ql*Q2/2,Q2xxx=-Qlx*C!2/2-C!l*C!2x/2, 

Qlxxxxx=-24*Qlx*qixx-8*qi*qixxx-q2x*q2xx*3/2-Q2*C!2xxx/2, 
Dq2=q2x , qixxx=-8*qi*qix-Q2*q2x/2 , ppp) ) : 

Pres:=coeff s(rrr, [L,Z] , 'tres') : 



printdif res : =proc (sss) 

local rr; 
for i from 1 to nops([tres]) do 
rr : =expand(Pres [i] ) : 
print ( ' nuinber= ' ) ; 
print (i) ; 

print (' spectral pareineters=' ) ; 
print (tres [i] ) ; 
print ( ' eleinent= ' ) ; 
print (rr) ; 

od; 
end: 



print ( ' spectral pareimeters ' , -tres [2] *162) ; 

EE : =expand (subs (q2xx=-qi*Q2/2 , Qlxx=-4*Ql"2-Q2"2/4 ,Pres [2] ) ) : 

EE:=expand(EE)/(81*2) : 

print (' energy is equal to res2'); 

print (EE) ; 

print ( ' spectral parameters ' , 3*tres [4] /4) ; 

K : =expaiid(subs (Q2xx=-Ql*Q2/2 , Qlxx=-4*Ql"2-Q2"2/4, Pres [4] ) ) : 
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K:=-4*expand(K)/3: 

print (' second integral is res4'); 

print (K) ; 

print ( ' spectral parameters ' , tres [7] ) ; 

res7:=expand(subs(C12xx=-C!l*C!2/2,Qlxx=-4*CJl"2-CJ2~2/4, 

C!lxxx=-8*Ql*Qlx-Q2*Q2x/2,Q2xxx=-Qlx*C!2/2-C!l*C!2x/2, 
Pres[7])): 
res7:=expand(res7) : 
print (' element res7='); 
print (res7) ; 

print ( ' spectral parameters ' , tres [3] ) ; 
res3:=expand(subs(q2xx=-Ql*Q2/2,qixx=-4*Ql'2-Q2'2/4, 

Qlxxx=-8*Ql*Qlx-Q2*Q2x/2,Q2xxx=-Qlx*Q2/2-Ql*Q2x/2, 
Pres[3])) : 
res3 : =expand(res3) : 
print (' element res3 ='); 
print (res3) ; 

print ( ' spectral parameters ' , tres [1] ) ; 

resl : =subs (Q2xx=-Ql*Q2/2 , Qlxx=-4*Ql~2-Q2~2/4 , 

Qlxxx=-8*Ql*Qlx-Q2*Q2x/2,Q2xxx=-Qlx*Q2/2-qi*C!2x/2, 
Pres[l]) : 
resl : =expand(resl) : 
print (' element resl='); 
print (resl) ; 

print ( ' spectral parameters ' , tres [5] ) ; 

res5 : =subs (Q2xx=-C!l*Q2/2 , Qlxx=-4*Ql~2-Q2"2/4 , 

qixxx=-8*Ql*Qlx-C12*Q2x/2,C12xxx=-Qlx*C12/2-Cll*Q2x/2, 
Pres[5]) : 
res5 : =expand(res5) : 
print (' element res5='); 
print (res5) ; 
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print ( ' spectral parameters ' , tres [8] ) ; 
res8:=subs(Q2xx=-Ql*Q2/2,Qlxx=-4*Ql~2-Q2'2/4, 

Qlxxx=-8*Ql*Qlx-Q2*Q2x/2,Q2xxx=-Qlx*Q2/2-Ql*Q2x/2, 
Pres [8] ) : 
resB : =expand(res8) : 
print (' element res8='); 
print (res8) ; 

print ( ' spectral parameters ' , tres [6] ) ; 
res6:=subs(Q2xx=-qi*Q2/2,Qlxx=-4*Ql~2-Q2'2/4, 

Qlxxx=-8*Ql*Qlx-Q2*Q2x/2,Q2xxx=-Qlx*Q2/2-Ql*Q2x/2, 
Pres[6]) : 
res6 : =expand(res6) : 
print (' element res6='); 
print (res6) ; 



alias (q[l]=q[l] (x) ,q[2]=q[2] (x)) ; alias (wp=wp(x)) ; 
alias (u=u(x) ) ; 

Eq[l] :=diff (q[l] ,x$2)+q[2] -2/2+3*q[l] -2+a[0] *q[l] -a[l] /4: 
Eq [2] : =q [2] -3*dif f (q [2] , x$2)+q [2] "4*q [1] +q [2] -4*a [0] /4+a [4] /4 : 
wpl :=(4*wp~3-g2*wp-g3) ' (1/2) : 
wp2:=eval(diff (wpl.x)) : 

wp3 : =eval (expand ( subs (dif f (wp , x) =wpl , wp2) ) ) ; 

F:=-6; B:=0; a[0]=0; G:=-144; 

for i from 1 to 2 do 

SEq [i] : =eval (expamd (subs (q [1] =F*wp+B , 

q[2]=(G*wp-2+E)-(l/2),Eq[i]))) : 

SSEq[i] :=expand(subs(diff (wp,x)~2=(4*wp~3-g2*wp-g3) , 

dif f (wp , x$2) =wp3 , expand (SEq [i] ) ) ) : 
#print('SSEq(' ,i, ') : = ' ,SSEq[i]) ; 



24 



n[i] :=degree(SSEq[i] ,wp) : 

for j from to n[i] do 

Cowp[i, j] : =coef f (SSEq[i] ,wp, j) : 

#Cowpl [i , j] : =expand(subs(F=0,G=0,Cowp [i , j] ) ) ; 

od: 

#FFl[i] :=seq(Cowp[i,il]=0,il=0. .n[i]) ; 
FF2[i] :=seq(Cowp[i,i3] ,i3=0. .ii[i]) ; 
od: 

for il from 1 to 2 do 
for jl from to n[il] do 

#print('Cowp(',il,',',jl, '):=', Cowp [il , j 1] ) ; 

printCCowpC ,il, ' , ' ,jl, ') : = ' ,subs(a[0]=0,E=36*g2, Cowp [il , j 1] ) ) ; 

od; 

od; 

#for i2 from 1 to 2 do 
#for j2 from to n[i2] do 

#print('Cowpl(' ,i2, ' , ' ,j2, ') : = ' ,Cowpl [i2, j2] ) ; 

#od: 

#od: 

#Fs:={FFl[l] ,FF1[2]>; 

#sols:=solve(Fs,{F,A,B,G,D,E,Cl~2,C2~2,gama,beta}) ; 
with(grobner) ; 
Fg:=[FF2[l] ,FF2[2]] ; 

#Fgl:=subs(eta-(-l)=etal,Cl"2=CCl,C2~2=CC2,Fg) ; 
#gsolve(Fg) ; 

#gbasis(Fgl , [A,b,F,G,D,E,CCl,CC2,r , s, sigma,gama,beta,eta,etal ,g2,g3] ,plex) ; 

#gbasis(Fgl, [A,D,Cl~2,C2~2,gama,beta] ,tdeg); 
#gsolve(Fgl, [A,D,Cl~2,C2"2,gama,beta] ) ; 

#Al:=expand(subs(A=(-6*s-sigma*D)/eta,B=A'2/F-F*g2/4,G=-eta*F/sigma, 

# E=D-2/G-G*g2/4,Cowp[l,3])) ; 

#D1 :=expand (subs (A=(-6*s-sigma*D) /eta, B=A'2/F-F*g2/4,G=-eta*F/sigma, 

# E=D-2/G-G*g2/4,Cowp[2,3])); 
#C12:=solve(Cowpl[l,0] ,Cl-2) ; 

#Gamal : =solve (Cowpl [1 , 1] .gaima) ; 
#Gama2 : =solve (Cowpl [1 , 2] ,gama) ; 
#AA : =solve (Cowpl [1 , 3] , A) ; 
#C22:=solve (Cowpl [2,0] ,C2-2) ; 
#Bel :=solve (Cowpl [2,1] ,beta) ; 
#Be2:=solve (Cowpl [2,2] ,beta) ; 
#DD:=solve (Cowpl [2,3] ,D) ; 
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